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In this report, we find the Mbh estimated from the formalism of Wang et al. (2009)[1] are more consistent with those from the 
Mbh-o", relation than those from previous single-epoch mass estimators, using a large sample of AGNs. Furthermore, we examine 
the differences between the line widths of H/3 and Mg II in detail by comparing their line profiles. The flux around the line core and 
that in the wing of both H/3 and Mg II show an opposite variation tendency, which indicates the BLR is multi-componential. The 
contribution of the wing makes the FWHM deviate from o"/,,,,,, and thus bias the Mbh estimated from previous single-epoch mass 
estimators. Thus the correction on the formalism suggested by Wang et al. (2009)[1] is crucial to Mbh estimation. 
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1 Introduction 

Accretion onto super-massive black holes (SMBHs) is gener- 
ally considered as the energy engine of active galactic nuclei 
(AGNs). The determination of the mass of SMBH (Mbh) is 
crucial to the understanding of most physical processes asso- 
ciated with SMBH and the cosmological evolution of black 
holes. The Mbh of type I AGNs are usually measured us- 
ing the virial theorem, MBH=fRBL«V~/G, if the size of broad 
line region (^blr) and the virial velocity (V) of clouds in the 
BLR are known, where f is a factor of order unity depend- 
ing on the geometry and kinematics of the BLR. Rblr can be 
estimated using the reverberation mapping (RM) method[2], 
which monitors the variability of continuum and emission 
lines. V can be estimated from the widths of emission lines. 
Conversely, Mbh can also be estimated using the tight cor- 
relation between Mbh and the stellar velocity dispersion of 
the galactic bulge (Mbh-o", relation)[3,4,5]. However, both 
of these methods cannot be used for large samples of AGNs, 



*Con'esponding author (email: wangjg@ynao.ac.cn) 



because the RM method is time-consuming and the measure- 
ments of 0", are limited by the spectral and saptial resolution 
of telescopes. 

For large samples of AGNs, Mbh can be estimated by 
combining Rblr, which is estimated using the important re- 
lationship between Rblr and the monochromatic continuum 
luminosity (R-L relation)[6,7,8], and the FWHM of emis- 
sion lines. The single-epoch mass estimators have been 
studied for various broad lines, such as HyS [1,9], Ha [10], 
Mgll /12800 [1,11,12] and CiV /11549[13]. If both H^ and 
Mgll FWHMs are good tracers of the virial velocity and 
can be used to estimated the Mbh, they should give the 
same Mbh values. Some researchers found they are con- 
sistent with each other [11,14,15,16], while others came to 
an opposite conclusion [1,17,18]. Wang et al. (2009) found 
that Mg II FWHM is systematically smaller than Hy8 FWHM, 
and that the relationships between WjS and Mg II FWHM and 
cr/,„(., which is the best virial velocity tracer measured on the 
variable part of the spectrum[19], deviate from the 1:1 re- 
lationship. The dependance of Mbh on FWHM should be 
Mbro'^FWHM'', where y is smaller than 2 for both Hy6 and 
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Mg II. If this is the case, most previous single-epoch mass es- 
timators (Mbh<^FWHM^) would introduce systematic biases 
in Mbh estimations[l, 20,21, 22,23] and result in many artifi- 
cial conclusions, as discussed by Rafiee and Hall (2011)[21] 
and Croom (20 11) [22]. Thus, further testing the validity of 
the formalism of Wang et al. (2009)[1] is critical for elimi- 
nating such biases in Mbh estimations and many other related 
relationships in AGNs. The Mbh estimated from the Hy6 and 
Mg II formaUsms of Wang et al. (2009) are more consistent 
with those from RM measurements than those from previous 
single-epoch mass estimators and are consistent with each 
other for a large sample culled from Sloan Digital Sky Survey 
(SDSS) Data Release 5 (DR5). However, one remaining issue 
is whether the new Mbh estimates are consistent with those 
derived from the Mbh-o", relation, which should be tested us- 
ing a large sample. 

Moreover, the reasons for the systematic deviations be- 
tween Hj6 and Mg II FWHM and cru„e are unclear now. The 
profile of an emission line is determined by the structure and 
kinematics of the BLR, which are complex. It is possible that 
the broad lines in most AGNs are generated in multi-regions, 
including the gravitationally-bound BLR, outflows [24] and 
the surface of accretion disk(Wang et al. 2005 [25]; Wu et 
al. 2008[26]). Different measurements of line width, such 
as (Tiine and FWHM, would represent different information 
about the structure and/or kinematics of the BLR. Special at- 
tention must be noted when using in the estimation of Mbh- 
The study of the structure and kinematics of the BLR would 
be helpful to understand why FWHM deviates from cr;,„g and 
important for the Mbh estimation of AGNs. In this report, 
we examine whether there are systematic biases between the 
Mbh estimated from the single-epoch mass estimators and 
those from the Mbh-o"« relation. We also compare the pro- 
files of Hy6 and Mg II in order to understand their differences 
and why their FWHM deviates from cr/,,,;.. 

2 The Bias in the Mbh Estimates 

We first verify the consistency between the Mbh estimated 
from the single-epoch mass estimators and those from the 
Mbh-ct. relation (Gultekin et al. 2009)[27]. We select 8470 
AGNs with z < 0.35 from SDSS DR4. The spectrum is cor- 
rected for the Galactic extinction using the extinction map 
of Schlegel et al. (1998)[28] and the reddening curve of 
Fitzpatrick (1999) [29]. The fitting method is described in 
Dong et al. (2008)[30] and described below. In the wave- 
length range 4030-7500A, we fit simultaneously the feature- 
less continuum and the Fell multiplets and other emission 
fines. Each of the [Olll] AA 4959,5007 doublets is mod- 
eled with two Gaussians, one for the line core and the other 
for the possible blue wing. The narrow components of Ha 
and HyS are fitted with similar profile to the line core of 
[Olll] /15007 and the broad components of them are fitted 
with 1-4 gaussians. Mbh can be estimated using the width 
of the line core of [Olll] /15007 as substitute for cr, [31]. Be- 



cause H/? is weak for many objects, we estimate H/3 FWHM 
from He FWHM[10] and then estimate the Mbh using the 
single-epoch mass estimators. We find that the differences 
between the Mbh estimated from the single-epoch mass es- 
timators and those from the Mbh-o"* relation are correlated 
with Ho- FWHM (Figure 1), if the formalisms from Greene 
and Ho (2005; hereafter GH05)[10] or Vestergaard and Pe- 
terson (2006; hereafter VP06)[9] are used. The correla- 
tion would decrease largely, if the formalism of Wang et al. 
(2009)[1] is adopted. The relationship between the Mbh dif- 
ferences and FWHM(Hq') is somewhat linear in log-log space 
and can be expressed as log ^ = klog ™^i""' +b. The 
(k,b) for the formaHsms of GH05, VP06 and Wang et al. 
(2009) given by the regression method of Kefly (2007) [32] 
are (2.05+0.04, -7.51+0.13), (1.93+0.03, -6.73+0.09) and 
(0.86+0.03, -2.99+0.11), respectively. All the intrinsic scat- 
ters of these relations are around 0.02 dex. This indicates 
that the Mbh estimated from the formalism of Wang et al. 
(2009)[1] are less biased than those from previous single- 
epoch mass estimators (Mbh^^FWHM^). We attempt to es- 
timate the Mbh using the Mbh-c, relation from other authors 
(Xiao et al. 2011 [33]) and find the Mbh estimated from the 
formalism of Wang et al. (2009)[1] are still less biased than 
those from previous single-epoch mass estimators. 




3.0 3.5 4.0 

log FWHM(Ha)/km s"' 



Figure 1 Correlations between FWHM(Hq') and the differences of Mbh 
estimate from single-epoch mass estimators and those from Mbh-o"« I'ela- 
tion[27] for the sample from SDSS DR4. The crosses are the median values 
and standard deviations of the Mbh differences and FWHM in each bin of 
FWHM. The solid lines show the best-fit relations. 



3 Profiles of and IMg ii 

The comparison above shows that the method of Wang et al. 
(2009)[1] is capable of correcting the systematic biases in the 
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Mbh estimations over a large redshift interval. This indicates 
indirectly that Hy6 and Mg II FWHMs are deviating from cr;,„p 
systematically. The systematic deviation may be caused by 
the complex structure and kinematics of the BLR. We com- 
pare the profiles of Hy6 and Mg II using the sample from Wang 
et al. (2009)[1], which was selected from SDSS DR5. The 
sample includes 495 AGNs with high signal-to-noise ratio 
(S/N > 20) in both the H/? (4600-5 100 A) and the Mg II (2700- 
2900 A) regions, which makes it suitable for the comparison. 
The spectrum is corrected for the Galactic extinction using 
the extinction map of Schlegel et al. (1998)[28] and the red- 
dening curve of Fitzpatrick ( 1 999) [29] . The redshifts of these 
quasars are from Hewett and Wild (20 10) [34], which were 
derived by cross-correlating observed spectra with a carefully 
constructed template. The dependence of emission line shift 
on luminosity and redshift are corrected and the systematic 
errors of redshifts are reduced to the level of 30 km/s, which 
are important to our investigation. We perform the contin- 
uum and emission-line fitting using an Interactive Data Lan- 
guage (IDL) code based on MPFIT [35], which performs 
minimization by the Levenberg-Marquardt technique. 
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Figure 2 BAI distributions of H/J (solid line) and Mg II (dotted line). 

The spectrum is fitted in two wavelength range: H/J range 
(4200-5600A) and Mg II range (2200-3500 A). For the Hj6 
range, the fitting method is similar to that described above. 
For the Mg II range, the method is described in Wang et al. 
(2009)[1]. The featureless continuum and Fell multiplets 
were modeled simultaneously. The broad component of each 
of the Mg llAA 2796,2803 doublets is modeled with a Gauss- 
Hermite series profile and the narrow component of each of 
the doublets is modeled by a Gaussian. 

Usually, the shift and asymmetry of lines are studied sep- 
arately, while they may be caused by the same process [24]. 
The blueshift and asymmetry index (BAI), which is defined 
as the flux ratio of the blue part to the total profile, measures 
their combined effects [24]. For H/? and Mg II, the blue part is 
the part at wavelength short than 4862.68 and 2800.26 A, re- 
spectively. The distributions of BAI are showed in Figure 2. 
The median value of the BAI distribution of Mg II is around 



0.5, while that of HyS is smaller than 0.5. Because both HyS 
and Mg II show no evidence of shift (see Figure 3), the BAI is 
primarily caused by the line asymmetry. This indicates that 
Mg II profile is quite symmetrical in that there are more flux 
in the red part of Hy6 than that in the blue part, which are con- 
sistent with the conclusion if the shift and asymmetry of lines 
are measured separately. 




Velocity (kn^/s) 

Figure 3 Composite profiles of H/3 and Mgll of four sub-sample 
divided by their Yifi FWHM, as well as their difference in veloc- 
ity space. First two panels are H/J and Mgll profiles. All these 
flux are normalized at the emission-line-free window 3030 - 3090 A 
and the continuum and Fell multiplets were subtracted. The last 
panel shows the line ratios of H/3 and Mgll. Black: FWHM(H/?) < 
3000 km/s; Red: 3000 km/s < FWHM(H/3) < 4000 km/s: Green: 
4000 km/s < FWHM(H/3) < 5500 km/s; Cyan: FWHM(H/3) > 5500 km/s. 

A direct comparison between the profiles of H/3 and Mg II 
is showed in Figure 3. The spectra are normalized at the 
emission-line-free window 3030 - 3090 A and the feature- 
less continuum and Fe II multiplets are subtracted. The sam- 
ple is divided into four sub-samples according to YifS FWHM. 
The composite HyS and Mg II profiles of each sub-sample, as 
well as their line ratio, are showed in Figure 3. The peaks of 
both HyS and Mg II do not show evident shift. The flux in the 
wings increases with the increase of FWHM, while the flux 
around the line core decreases with the increase of FWHM. 
The change of H/? is more rapid than that of Mg II. This indi- 
cates that Hj6 and Mg II are not cospatial in BLR and the BLR 
in AGNs is multi-componential. At least two emitting regions 
are needed: an intermediate line region (ILR) producing the 
line core and a very broad line region (VBLR) producing the 
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line wings [36]. The emission in ILR makes larger contribu- 
tion to the Mg II lines, while the emission in the VBLR makes 
larger contribution to the Balmer lines[36]. 



4 Discussion 

Different structures of the BLR have been proposed to ex- 
plain the profiles of emission lines. These models includes: 
a rotating accretion disk, binary black holes, bipolar outflow 
and anisotropically illuminated spherical BLR (see Eracleous 
and Halpern 2003 and reference therein) [37], as well as the 
gravitationally-bound BLR-Houtflow model of Wang et al. 
(20 11) [24]. The gravitationally-bound BLR+outflow model 
has succeeded in explaining the profiles of the high ioniza- 
tion C IV line [24], but is not suitable to explain the profile of 
low ionization H/? line. This is because H/? shows a system- 
atically small BAI (<0.5) opposed to the expectation of the 
model (BAI>0.5). Eracleous and Halpern (2003) found the 
accretion disk emission could explain the double-peak pro- 
file and other spectroscopic properties of AGNs presenting 
the double-peaked Balmer lines, while other structures are 
unsatisfactory[37]. They attempted to explain the profiles of 
HjS and Mg II using the accretion disk emission, but the model 
predicts lower Mg II flux than the observed flux. One of the 
possible reasons is that the BLR is two-componential. The 
contribution of the ILR to Mg II is critical but is not included 
in their model. 

As showed in Figure 3, the contribution of the VBLR to 
Yip and Mg II flux becomes more important with the increase 
of Hj6 FWHM. However, the contribution of the VBLR to 
(Tiine may be small, because clouds in the VBLR might be 
optically thin to the ionization continuum[38]. The emission 
from the VBLR does not vary with the variability of contin- 
uum and contributes little to the variable part of the spec- 
trum. This may be the reason of the systematic deviations 
between H/? and Mg II FWHM and crune- Morevoer, the frac- 
tion of the contribution of the VBLR to Mg II is much smaller 
than that to Hy8, which makes the Mg II FWHM systemati- 
cally smaller than H/3 FWHM. The contribution of the VBLR 
makes FWHM deviate from cr;,-,,^ systematically and bias the 
Mbh estimation from previous single-epoch mass estimators 
(Mbh ocFWHM^). When estimating the Mbh using FWHM 
as the tracer of the virial velocity, it is crucial to correct the 
biases by using the fitted index of the Mbho^FWHM^ relation, 
rather than the assumed y = 2, as suggested by Wang et al. 
(2009). 
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